Cosmological perturbations of Quasi-Dilaton non-linear massive gravity 
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We study the cosmological perturbations of the recently proposed extension of non-linear massive 
gravity with a scalar field. The added scalar field ensures a new symmetry on the field space of the 
theory. The theory has the property of having a flat dS solution, in contrast to the standard dRGT 
massive gravity. The tensor part is the same as that of the standard dRGT massive gravity and 
shows the gravitational wave with time dependent mass which modified the dispersion relation. We 
obtain the curvature perturbation of the model on superhorizon scales for a specific choice of ui — 
and find that there is a possibility that the curvature perturbation is damped over such scales. 

PACS numbers: 



I. INTRODUCTION 



Gravity, as the oldest force known to man has also had the longest history of trials and tribulations along the road 
£N) • to discovering its nature. Over the course of its development, it has been witnessing a myriad of attempts to unlock 
its notoriously difficult and mysterious behavior from the largest to the smallest of distances. The challenge has been 
truly spectacular. Since the first formulation of the gravitational field by Newton and centuries later by Einstein in 
the form of the theory of general relativity (GR), the scenery is still cluttered with debris left from various attempts to 
understand its hard to grasp nature. Even today, the challenge is as fresh and as interesting as ever. Not surprisingly, 
building on GR, the last couple of decades have been particularly rich in new ideas and approaches which attempt to 
formulate the gravitational field in such a way as to pave the way to a formulation of the theory which would explain 
i— such recently observed phenomenon as the accelerated expansion of the universe, galaxy rotation curves and even 
the birth of the universe. One such attempt has been surfacing over the past few years in the form of what is now 
^ ', known as massive gravity which, as the name suggests, is a theory with a massive graviton as the building block of 
$H ' the gravitational field. 

The notion of a massive graviton has been a tempting and challenging premise in theoretical physics. One of the 
main motivations of having a massive graviton is that gravity could become weak at large distances, thus mimicking 
-— I ' the effect of accelerated expanding universe. The problem is not as easy as it seems. The first attempt to build a 
^ ■ theory for a massive spin-2 field was back in 1939 when Fierz and Pauli (FP) developed a linear theory for a massive 
graviton It took thirty years for physicists to find out that the theory does not reduce to the standard GR 
when one takes the limit m — >• 2}. The problem was soon addressed by Vainstein [U who proposed that adding 
non-linearities to the action can cure the problem and screen the effect of helicity-0 component of the massive graviton 
at solar system scales. The simplest possible non-linearity one can add to the FP action is by replacing the linear 
kinetic term for the helicity-2 field with the fully non-linear, and still ghost free, Einstein-Hilbert action. However, 
the resulting action have proven to have a ghost instability which was discovered by Boulware and Deser [4]. The 
problem arises because the lapse function is no longer a Lagrange multiplier. This new problem can then be solved if 
one appropriately adds interaction terms to the Lagrangian and again makes the lapse function a Lagrange multiplier 
order by order in non-linearities 0, 0] . 

In this paper, after a brief review of the theory we study the cosmological perturbations around an accelerating 
solution and obtain the background equations. The second order Lagrangian can then be obtained using the perturbed 
■ FRW metric. It is an immediate observation from the form of the second order Lagrangian that the tensor, vector 
and scalar modes do not couple. The tensor mode shows a massive gravitational wave with a time-dependent mass 
parameter which we shall obtain in section IIV Al We show in section IIVBI that the vector mode vanishes when one 
integrates out the non-dynamical vector field and then obtain the scalar mode in section IIV CI Two of the scalar 
modes are non-dynamical and can be integrated out immediately with three degrees of freedom remaining in the 
scalar sector. We show that in the superhorizon limit where one of the degrees of freedom does not play a role, the 
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curvature perturbation can be obtained analytically, and there is a vast region in the parameter space over which the 
curvature perturbation will damp on the superhorizon scales. 



II. A BRIEF REVIEW 

Recently, a two parameter theory of massive gravity has been proposed 0] • With the aid of a new re-summation 
of the non-linear interaction [8| de Rham, Gabadadze and Tolley (dRGT) have constructed a theory which is free of 
ghost instabilities in the decoupling limit. The theory has been proven to be ghost-free in the full non-linear theory 
and hence a reliable effective field theory of a massive spin-2 field [9] . The action for dRGT non-linear massive gravity 
can be written as 



where the non-linear interactions are collected in U 

U{K) =U 2 +a 3 U 3 + a 4 U 4l (2) 



which consists of polynomials of various traces of the matrix tC£(g, <fi a ) = 5% — \J g^ a ]av where the fiducial metric is 
defined as f av = d a 4> a d„(j> h r]ab and <j) a are the Stuckelberg fields responsible for the breaking of general covariance and 
the rectangular brackets denote traces, [K] = Tr(/C) = /C^, 

U 2 = [JCf - [K 2 ] , (3a) 
W 3 = [/C] 3 -3[/C][/C 2 ]+2[/C 3 ], (3b) 
Ui = [JCf - 6[/C 2 ] [JC] 2 + 8[/C 3 ] [K] + 3[/C 2 ] 2 - 6[/C 4 ] . (3c) 

The first term concides with the Fierz-Pauli mass term at the linear level, and the last two terms are non-linear 
interactions which ensure that the theory has no ghost. 

One of the interesting property of the theory is that if one assumes the Stuckelberg fields to be in the unitary gauge 
defined as <fi a = <^x M , the theory does not have a non-trivial flat FRW solution [10(. However, the theory has an open 
FRW solution with an additional consideration in that one must transform the field space to the open slicing of the 
Minkowski metric [ll|. The cosmological perturbations around such a solution was also considered in [l2| where the 
authors found that the scalar, vector and tensor modes actually decouple and as a result the vector mode does not 
play any role in the theory. The tensor mode describes a massive gravitational wave with a time-dependent mass. 
Cosmological evidences of dRGT theory is also considered in 13]. For a review on the theoretical aspect of the theory 
see 

One can also let the non-dynamical metric of the theory to have a kinetic term and hence construct a bimetric theory 
15| which has been proven to be ghost free [l6j]. The cosmological aspects of such bimetric theory is investigated in 
17j j. One may continue the procedure to build a multimetric gravity theory with dRGT non- linear mass terms. In 18], 
the authors show that the theory is ghost free in the metric formulation only when the interactions between gravitons 
are not cyclic. However, in [191 ] the authors show that in the vielbein formulation of the theory any interactions 
are allowed. Also in [2(| the authors show that the problem of having no flat FRW solution persists in multimetric 
theories. In fact if one of the metrics is assumed to be flat, all of the other metrics will become flat. 

One of the solutions to the problem of the non-existence of the flat FRW solution in the theory is by extending 
the theory in such a way that the graviton mass becomes a function of some scalar field ip [2lJ. The cosmological 
solutions and dynamical analysis of such theories are considered in [22| . 

Another way to extend the theory is to couple a scalar field to the mass Lagrangian such that the resulting new 
Lagrangian has an extra symmetry. In particular, one can couple the scalar field to the Lagrangian to achieve dilation 
invariance on the field space [23| 

cr a- M pl a, 4> a -> e a (f> a . (4) 
In the Einstein frame one can write the new action as 



— ^/ — n hi — 

M, 



R - —g^d^d v a + m 2 U{t) 



+ V zz g^m(g uu ,ip), (5) 
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where /C£J is now denned as 

K{g, <f> a ) =SH- e^' (6) 

Only the purely geometrical part of the above action is invariant under transformation (U), thus the acronym Quasi- 
Dilaton (QD) for the scalar field [23J. This theory has been proven to be free of ghost in the Minkowski background 
if uj > 6. The most interesting feature of this theory is that it admits a flat de Sitter solution even if the Stuckelberg 
fields are in the unitary gauge [23j]. The de Sitter solution is also stable in the decoupling limit for 

a 3 ^0, 0<a 4 <-^, 0<w<6. (7) 



III. THE BACKGROUND EQUATION 

Let us assume that the background metric is of the form 

ds 2 = -N{tfdt 2 + a(t) 2 (dx 2 + dy 2 + dz 2 ) , 
and the Stuckelberg fields take the form 

0° = /(*), t l = 5W. 



(8) 



Note that we will work in the unitary gauge. However, in order to obtain the Stuckelberg equation from the action 
we assume the above form for the Stuckelberg fields, and finally set f(t) = t. We also assume that the QD field only 
depends on t. By varying the action (TTJ) with respect to f(t), one obtains the constraint equation as 



2a/M pl 



-a± - a 3 | ae a/Mpl + I -a 4 + a 3 + 



(9) 



We are interested to the set of equations for which kx = 0. In this case one can then solve the above equation by 
using the ansatz 



e a/M pl = x 



in the equation which results in 



X 



3a3 + 8o!4 ± ^/9«3 — 16a!4 
8tt4 



X = 1. 



(10) 



The solution X = 1 is not acceptable because in this case the effective cosmological constant vanishes, and the 
consistency of the theory does not allow one to have a flat cosmological solution [23[ . Putting the above ansatz in the 
mass part of the Lagrangian, one can write the Friedman and Raychaudhuri equations as 



m 1 - 



M 2 vl 2N 2 



3M /t 



2§ + 3H 2 



uj a 
M 2 t 2N 2 = 3 9 " 



where we define H = dj (aN) and 



M t = 



m 

M a = — 

9 3 



16aJ 

2 



3aj (3(X ~ l)a 3 - 1) + 4 (1 - 4(X - l)a 3 ) a 4 



- 6 - X(-6 - 3r + X + 2rX) + 3{X - 1)(4 + (-2 + r(X - 3))X)a 3 + 12(X - l) 2 (rX - l)a 4 



(11) 
(12) 

(13) 
(14) 



and r = r(t) — a(t)/N{t). We will use these forms of the background equations in the second order Lagrangian. 
Substituting equation (fTTJ)) in (fTTj) leads to a constant Hubble parameter, showing a de Sitter solution 



H 2 = H 2 



'6m 2 (l - X) [2 + 4a 3 + 4a 4 + X(-l + (X - 5)a 3 + i(X - 2)a 4 )] 



1/2 



uj — 6 



(15) 
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From equation ([12")) we obtain that r should be a constant given by 

_ 12(1 + 2a 3 + 2a 4 )w + 4X 2 (1 + 6a 3 + 12a 4 )(3 +uj) - 3X 3 (a 3 + 4a 4 )(6 + w)- 3X(1 + 3a 3 + 4a 4 )(6 + 5w) 
X(u - 6)(3 + 9a 3 + 12a 4 + 3X 2 (a 3 + 4a 4 ) - 2X(1 + 6a 3 + 12a 4 )) ' 

(16) 

The scalar field equation is satisfied automatically by plugging in equations (|10p. (fl5|) and (|16p . 

IV. SECOND ORDER LAGRANGIAN 

In this section we study perturbations signified by 

ds 2 = -N(t) 2 [l + 2c/)(t,x,y,z)]dt 2 + 2N(t)a(t)/3 t (t,x,y, z)dtdx l + a(t) 2 [5 tj + h zj (t,x,y,z)]dx l dx j , (17) 

where 4>, and hij are the perturbation variables of the FRW metric. The perturbation of the Stuckelberg fields 
around unitary gauge is 

r = x a +7r a (t,x,y,z), (18) 
and the perturbation of the dilaton field around the background solution has the form 

a = <T {t) + ((t,x,y,z), (19) 
Now, we consider the infinitesimal coordinate transformation 

as" -> x» +^(t,x,y,z), (20) 
which leads to the change of the perturbed quantities as 

«^0-l9 t (jV£ ), (21a) 

Pi Pi + -did - ^M, (21b) 

hij ->■ /»„ - ~ fyfi ~ 2NHt%j , (21c) 
7r° ->• 7r a - (21d) 
C^C-^ , (21e) 

where prime denotes the time derivation. 

Using the perturbations of the Stuckelberg fields one can construct the gauge invariant quantities using the perturbed 
Stuckelberg fields in the following manner 

(22a) 
(22b) 

NH7r°Sij, (22c) 

(22d) 

We may decompose the gauge invariant vector and tensor parts as 

B l = diP + S h (23a) 
Hij = 2^S tJ + didjE + ^(diFj + djFi) + (23b) 

where 

d'S, = = 9\F i; 

a i 7«=0 = 5 iJ 7y. (24) 





= 4>- 




B, 


= Pi- 


a N 


Hij 


— 


- di-Kj - djiTi 


Z 


= c- 
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We note that there are four degrees of gauge freedom because of the coordinate transformation, two of which represent 
the scalar part and the others relate to the vector part. One may fix the gauge freedom by the choice 

7T° = 0, 7T 1 = 0. (25) 

Note that with this gauge fixing all the gauge invariant perturbation variables become equal to the original one, 
making our calculations simpler. Also note that the above gauge fixing is similar to the use of the unitary gauge, and 
so the fiducial metric takes the form 

f^u = V- ( 26 ) 
The components of the fjf — j pv matrix in the unitary gauge are 

/o° = ^ (1 - 20 + 40 2 - A/3 1 ) + 0(e 3 ), (27a) 

fi=Jfo (A ~ P* h Ji ~ 2 ^) + °( £3 )> ( 27b ) 
fo = (P l ~ Pjh Si - 2#T) + 0(e 3 ), (27c) 

fj = ^ (<Sj - ti s - Ffa + h ik h kj ) + 0(e 3 ). (27d) 

To compute the components of the /C£f, we use the method presented in [12j to expand the square root in (pj). 
To zeroth order perturbation we find 

£ (0)0 = 1-^, £ (0)i = = K (0)0 4 , ^.= (1-^5^ (28) 



where A = e a °^ Mpl . The first and second orders are 



' N \ M pl J' ^ l ~ N{l + r) Pu ^ °~N(l+rf' ^ * a ^ M pl ' 



~ ( 2)o A f r(r + 2) , 3 2 1 J_ 2 \ 

* = n [wrw A " ^ + ^ C0 ~ 2a*j c J ' (30a) 

^ (2 \° - TTT^T ( + - -S-CjfO , (30b) 



AT(r + l)\r+l 2(r + iy J1 M pl 

where r — jSj , as that in the background. The traces of K for zero and first order perturbation are given by 

pC™] (0) =3(l-X) n + (l-rX) n , (31) 



[«•]<» - -Y,l - rX)-' -^C) + =X(1 - X)"- (* - JU) , 



(32) 
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where X = A/a. To second order we obtain 



[^2] (2) = r2 Xr3 XF0i + {ArX - 3)A>0 2 + - 2rX)Xr0C + t^t [ 3 ( 2 ^ ~ 1) + r ( 2rX - l )] X C 2 

T\ Mpi M pl 

+ -J-(l - 2X)X/< + (X - ^:.\7,'</, (f . (33b) 
[^ 3 ] (2) = TT-fa ~ 2r 3 X + r 4 X 2 )X/3 l ft - |(3 - 5X)(1 - X)Xh ij h ij - |(1 - rX)(3 - 5rX)Xr</> 2 
+ - rX)(l - 3rX)Xr<K + ^-(1 - - 3X)X/< 

° [(r + 3)-4X(r 2 + 3)+3X 2 (r 3 + 3)]XC 2 , (33c) 



[JC 4 ]( 2 ) = A ( r2 _ 3Xr 3 + 3X 2 r 4 - X 3 r 5 ) X(3 % (3 1 + |(2X - 1)(1 - X) 2 Xhijh ij + 6{2rX - 1)(1 - rX) 2 Xr</> 2 
+ -jj7-(l - 4rX)(l - rX) 2 Xr<K + (-r - 3 + 6X(r 2 + 3) - 9X 2 (r 3 + 3) + 4X 3 (r 4 + 3)) X( 2 
+ -?-(! -4X)(l-X) 2 Xh(. (33d) 



where 



r n = 

i=0 



£V. (34) 



From the above formulae, one can construct the mass term using equations ((3J)- The gauge invariant second order 
Lagrangian can then be written as 



= Mpi J d 4 xNa 3 (c + ~M f (-$ 2 + B% + Sft) + ^(2M f - M g ){H 2 - 27/,///" ! + . 



(35) 



where we have defined 

C = ^ ill ' ll ' J - ^ + \ - 2^) ~ ^a 9 ^ ~ 3H ^ 2 + i iiW - n '' ~ {dlBl?) 

+ + mz + ^ zv2z + ^ zd,e ' + ■ (36) 

and 

C m ass = M ){,/){'■■ + M 2 H 2 + M C Z 2 + (M hC H + M&*)Z + M $ 2 + M p B t B l + (37) 

The definition of M,'s are given in Appendix [5] We have also used equations (fTTj) and (fT2"]) to simplify the action. 
The above action ensures that the scalar, vector and tensor modes do not couple to each other. Therefore we study 
them separately. 

A. Tensor mode 

Keeping only 7^ in the action (1351) . we obtain 

1 _ .. 1 

(38) 



SilLor = Kl I dixNa3 
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Variation of the above action with respect to 7y leads to the equation of motion for tensor perturbations 

1 - NaV 2 -y ij + 2(3Af g - AM^Na 3 ^ = 0. ( 3!)'l 

Fourier transforming the equation and using the conformal time define as 



dt\N 



one can write equation (|39l) as 



fc 2 - 



a?7 = — at, 
a 



2a 2 (3M g - 4 Afi) 



7 = 0, 



(40) 



(41) 



where we drop the indices of 7y and define 7 = ^7. This equation shows that the graviton acquires a time-dependent 
mass in this background. This is in agreement with the result of [l2j with a different mass parameter. 



B. Vector mode 



We now study the vector mode of the action There are two vector modes S l and F % in the action. One can 

write the vector part of the action as 



Sfl tor = M pl I d^xNa 



' F,V 2 F i + -±-F\7 2 S l - -Ks i V 2 S i 



16V 2 



4Va 



4a 2 



\{2,M f + 2M p )S l S l + -(3M S - 4M 1 )F l \7 2 F i 
2 8 



Varying the action with respect to S l gives 



1 



4Va 



V 2 F, - ^V 2 S, + {3M f + 2M )Si = 0. 



By inserting V 2 ^ to the action one obtains 



Suitor = ^QMg-Uf!) J d A xNo?F^ 2 F\ 



,(2) Mpl_ 
vector g 

which gives F{ = and then 5j = 0. So the vector mode of the action vanishes. 



(42) 



(43) 



(44) 



C. Scalar mode 



In this section we study the scalar perturbations of action (|35[) . We are interested in the behavior of perturbed 
variables at the superhorizon scales. The scalar part of the second order Lagrangian can be written as 



12M 



pi 



N 



i (Mi + M 2 - hl f + hlg) N 2 a 2 {V 2 E) 2 + -^N 2 \M K a 2 Z - ^V 2 $ ) V 2 E 



+ ^V 2 V>V 2 (V 2 S) + i yMf + ^M h ^j V 2 a 2 $V 2 £ + Qa/j + M 2 + ±M f - ^M g ) N 2 a 2 ^ (3^ + V 2 E) 
+ -^Na 2 H<S>V 2 E + io (aNH$ + ^VV 2 ^ - \ aW2 ^j ^ ~ \ a2 ^ 2 + \ ( M f + ~M h <s> ) X 2 n 2 <Pi 

_ -N 2 $\J 2 ib --( 2H 2 - -AU + Mf] a 2 V 2 $ 2 - -N 2 aH$S7 2 B + —M C t,N 2 a 2 Z<$> 
6 8 V 3 V 6 12 4 

+ X -M hQ N 2 a 2 Z^ - ^V 2 ^V 2 ^ + -^M ( N 2 a 2 Z 2 - -L (a^ + ~Ai>) V 2 a 2 /3V 2 /3 



1 



12A/ 2 V 2 



a 2 Z 2 - aW$ - 3V- 



^V 2 E)Z + ia 2 o- 2 $ 2 + iv 2 ZV 2 Z + o-oVaZV 2 ^ 



(45) 
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One can see from equation above that /3 and are non-dynamical. Transforming back to Fourier space, one can find 
their equations of motion 



-2M 2 HN<S> + ljZ& + 2M 2 i/> 
M 2 {2M p + 3M f )aN ' 



and 



k NHZ — a (6Mf + AMp)Z )a -6M 2 NH[ -k z - 6M f - <±M P )ip - M 2 N[ k 2 a z HE 



(46) 



+ a 2 N [ -M f + M h(j> ) (k 2 E - 6<tp) + a A M ((j> N( + 2k 2 'Nip 



where we have defined 



A = a 2 M 2 N 2 (6M f + 4Mp) (3M f - 2M rj> + 6H 2 ) + 8k 2 M 2 N 2 H 2 - a 2 uj(6M f + 4Mp)o 



(47) 



(48) 



At this point we are interested to the behavior of the fields on the superhorizon scales where k 2 — > 0. After substitution 
of (3 and $ from equations (pm) and (|T7|) and defining the curvature perturbation as 



K = ip + —Z, 

the Lagrangian (|4"5"j) at the superhorizon scales takes the form 

/ 1 



^scah 



&M 2 pl a 2 



- 3 [2M - 3M f + (oj- 6)H 2 ] \2 
1 



-r 2 (X 5 a 2 + 2rujHla + r 2 {2M^ - 3M f + ojH$)) V> 2 



1 



- r 2 (rujH 2 + \ 5 a)ak^ + -r 2 \ 5 a 2 K 2 - -A 3 a 4 (ft - ip) 2 - rX 2 a 3 ^(TZ - i>) - 2r 2 \ x a 2 ^ 2 



Hnr 



^X 6 a 2 {TZ -ip) + ra[X^ + (M p/ Af 0C - ujH 2 )TZ] + 9r 2 \^M H + M^j 



dip 



-H r 
2 ° 



- X 6 a(TZ -^)+ujr [(u - 6)H 2 + 2M 4> + 2M h( p] ij; 



a 2 tp 



(49) 



(50) 



where we define 



/3 3 \ 45/9 9 

Ai = \ -M } - -M g + Ah + 3M 2 J (w - 6) H 2 - —M 2 - f 3M X + -M H + 9M 2 - -M, 



, - ) M f 



6M 2 - -M g + 2M X M, 



-Mr 



A 2 = -u(u - 6)H 2 + y{M p iM K + M$ + M h4> )uj - 6M pt M h4 )H 2 
- 3M p i UM hC + ^M h<p )M f + ^Mh^M^ - ^M M K J , 



A 3 = -3ujH* + ((M + M 2 t M c + M pl M c<f> )uj - 6M 2 M c ) H 2 - 3M 2 ; Qm 2 + M c M f - ^M m c 



A 4 = -u{u - 4)ff 2 + (M w + M^u - M pl M <rt> , 



A, = 



1 



1 



H " ~ i KU + 2 Mf I W ' 



A 6 = -6A 5 + wMpiMtf. 
For uj = 0, the field equations at the superhorizon scales are simplified to 

3r 2 H MM c< pi> - a 2 \ 3 {TZ - ip) - 3raA 2 i/> = 0, 



(51) 

(52) 

(53) 

(54) 

(55) 
(56) 

(57) 
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^HqMM^TZ - \x 3 a 2 (TZ -ip) + 2ra(2\ 4 H% - \ 2 )ip 



r 3 (2M - 3M f )(r^ + 2aH ip) + a 

+ ra(4H MM Cq > + \ 2 )K + r 2 (9H 2 (3M f + 2M h<t> ) + 4Ai) V 



= 0, (58) 



which can be analytically solved for 1Z and ip. Substituting 1Z from (|57p into the second equation and solving the 
resulting equation for ij) results in 

/ \/ 9a2 h' 2 -32AB -^9A 2 h' 2 -32A \ 

ij) = ti I Cit 5 ™5 + C 2 t ^ , (59) 



where C\ and C 2 are integration constants and we have defined 

2 1 

3 * c 6 



A = M f M c - -M^M C + -M 2 ^ (60) 



J / . , ., . , . .. . r .. ■ ■ - • r .. - • r » ' 1 - • ^ t t2 , ^ / o /I » r2 7i f i nn7i/f2 I 



B =- i^Ah + 24M h - 6M h4> - 3M g - 3M f + 3(Af C0 - M hi )M K \ H£ + — I 2AM^M^ + 90MfM ( 
+ 24M^Af c - 24:M h4> M hC M C4> - (8Mj - 3M 9 + 24Af /l )(4M M c - Af c 2 ) 

+ 6M/ ( - 6M 2 ( + (8Mi - 3M g + 24M h + 12M, i0 - 4Af ) M c - 6M K M C4> + M 2 ^) ^ , (61) 



(^/9A 2 ff2-32AB -^/ 9 A 2 H{- -32A \ 

C 3 i + c* 4 t , (62) 

where C3 and C4 are some functions of C\, C 2 and Mj. To have a decaying mode for the curvature perturbation, one 
should keep the positive power of time only and this leads to 



a3 < & 04 < or —a 2 < 0:4 < ~^ a \ ff° r — 1 < 03 < — - we have 0:4 7^ — (1 + 3a3) ) , (63a) 
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a 3 > & a 4 < —a 2 3 , a 4 ^ 0. (63b) 

V. CONCLUSIONS AND FINAL REMARKS 

In this paper we studied the cosmological perturbations of the Quasi-Dilaton massive gravity. This theory is the 
extension of the non-linear massive gravity theory recently proposed by de Rham, Gabadadze and Tolley through 
a scalar field. The scalar field is coupled to the mass term in such a way that the field space of the theory admits 
a dilatation invariance. This new symmetry of the theory enables us to obtain flat FRW solutions. If considered 
without matter, the theory predicts an accelerating solution which is the effect of the graviton mass. The stability of 
this solutions is considered in [23| where the authors find that the u> parameter has to have a positive value less than 
6, and the range for the parameter a 4 has to be less than a|/8. 

Our analysis of cosmological perturbations shows that the vector mode vanishes when one uses the gauge invariant 
variables. The tensor mode has a different behavior to that of the standard GR but similar to the gravitational 
waves obtained in the dRGT massive gravity theory. The gravitational waves in this theory have a non-vanishing 
time-dependent mass which modifies the dispersion relation of the gravitational waves. 

In order to find the scalar spectrum of the theory, one can use the gauge invariant variables and then integrate out 
two of the non-dynamical variables included in the metric perturbation. The equations of motion of the remaining 
three scalar perturbations can then be obtained by varying the resulting action. As a matter of fact, the procedure 
is so difficult that one cannot solve the equations analytically. However, as long as we are interested in studying the 
behavior of the theory at the superhorizon scales, we can study the lagrangian over that scales. One of the scalar 
perturbations does not play any role in the superhorizon scales because it always comes with a wave number. The 
resulting superhorizon Lagrangian can then be varied with respect to ip and the curvature perturbation 1Z. The 
equations can be solved analytically if one assumes that the quasi-dilaton field has no kinetic term. We can then 
obtain conditions for which the curvature perturbation is damped over the superhorizon scales, i. e. equation (1631) . 
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However, if one considers the range for a 4 where the solution is stable (equation ([7}), one may reduce the allowed 
parameter space to be 



a 3 > and < a 4 < -p, 



(64) 



for uj = 0. One may note that the above range for the parameter space can be different if one considers a non-zero u> 
parameter. 

Finally, the QD massive gravity theory has a potential of producing a reasonable inflationary scenario, if one adds 
a potential to the action or change the graviton mass to be a function of the quasi-dilaton field. This can then break 
the dilatation invariance, which we will study in future works. 
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Appendix A: Constants of the second order mass Lagrangian 



Mi = 



128a| 



9raf(-l + 3(X - l)a 3 ) + 4 (l + 3r - 3(-3 + 4r(X - 1) + X)a 3 - 18(X - l)a§) a 4 



-16(8 + 9a 3 - X(7 + 6a 3 ) + 3r(-2 + X - 4a 3 + ZXa 3 ))a\ - 192(r - 1)(X - l)a\ 



(Al) 



M, 



32a 4 



r(l + 6a 3 )(3(A - l)a 3 - 1) + 4(6 - 4r - 3A - rX + 3(4 - 3A + r(4A - 5))a 3 a 4 + 48(r -1)(X- l)a 4 

(A2) 



M, 



C 32a 2 



3m 



9(7r - 3)a^(3(A - l)a 3 - 1) + 4( - 4 + 12r + 3(3 + r + 7X - 19rA)a 3 + 18 (9r - 5)(X - l)a^)a 4 



16(2 + X(U - 33a 3 ) + 39a 3 + r(6 - 51a 3 + 15X(3a 3 - 2)))a\ + 384(r - 1)(X - l)a\ 



(A3) 



3m 2 

M = — T [-3aJ + 9(A - l)a| + 4a 4 - 16(A - l)a 3 a 4 ] , 



(A4) 



M, 



hC 32a 2 



27ra 2 (l - 3(A - 1)0*3) - 4 (2 + 4r - 3(-6 + 3r + 2A + 7rX)a 3 + 18(-2 + 5r)(A - l)a|) a 4 



16(8 + X(2 - 24a 3 ) + 30a 3 + r(-2 - 39a 3 + A(-14 + 33a 3 )))a^ - 384(r - 1)(X - l)a\ 



(A5) 



M, 



in 



2 r 



C0 ~16a 2 



27(2(r - l)r - l)a 3 (3(A - l)a 3 - 1) +4 - 6 + 8(r - l)r + 3(9(A - 1) + 2r(2 + r + 7 A - 7rA))a 3 



(A6) 



+ 54r(2r - 3)(A - l)a\ a 4 + 16(4A - 6 + 8r(r -1 + 3X- 2rX) + 9(1 + 2(r - 2)r){X - l)a 3 )a\ 



Mm, = 



8a 4 



(1 + 6a 3 )(3(A - l)a 3 - 1) + 4(2 - 3a 3 + A(-4 + 3a 3 ))a 4 



(A7) 
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M a = 



32ria 4 



9ai n + 6rr 2 - 2(3 + 4r)r 3 + (8 + 3r)r 4 - 3r 5 3(X - l)a 3 - 1 



+ 4l 8rr 2 - (8 + 9r)r 3 + 3(3 + r)r 4 - 3r 5 + 3 



- 4{X - l)n - 2r(3 + 7X)r 2 + 6r 3 + 13rr 3 + 14Xr 3 



+ 17rA> 3 - 13r 4 - 6rr 4 - 17A> 4 - 6rXr 4 + 6(1 + A> 5 I a 3 + 6{X - 1) [lOrr 2 - 5(2 + 3r)r 3 



+ 3(5 + 2r)r 4 -6r 5 ]a^ 



|a 4 + 16 jrr 2 [ - 4 - 25AT + (75AT - 78)a 3 ] + r 3 [4 + 12r + 25X + 33rX 



+ 3(26 + 43r - 25X - 42rAT)tv 3 ] + 6r 5 [l + 2X- 9(X - l)a 3 ] 

+ 3r 4 ^ - 4 - 2r - 11X - 4rX +[-43 + l8r(X - 1) + 42X] a-^j \a\ 

+ 192 [r(5X - 4)r 2 + (4 + 7r - 5X - 9rX)r 3 + (-7 - 3r + 9X + 4rX)r± + (3 - 4X)r 5 ] a| 



(A8) 
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